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INTRODUCCION: 


Los seres organicos tienden a incrementarse con rapidez, motivo por el que de manera 
inevitable siempre hay una batalla por la existencia ya sea un individuo con otro de su misma 
especie o con los de alguna especie diferente o bien con las condiciones fisicas de la vida. La 
cantidad de alimento para cada especie da el limite extremo al que se puede aumentar cada 
una; pero muy a menudo no solo la obtencion de alimento si no servir de presa a otros animales 
es lo que determina el numero promedio de una especie.(Darwin,1882) 

Los modelos de ecosistemas resultan interesantes y utiles, tanto desde el punto de vista 
teorico como practico. La mayor parte de los estudios relativos a este asunto se centran en el 
desarrollo de herramientas que permitan predecir la evolucion futura de los ecosistemas 
sometidos a ciertas condiciones, con el fin de introducir tecnicas de control en los mismos 
(Volterra, 1931), (Lotka, 1952), (Dempster, 1975), (Braum, 1983). Tal es el caso del control 
biologico de plagas, es decir, el control de una especie que crece descontroladamente en un 
entorno o ecosistema para evitar que provoque danos colaterales. Elio como alternativa a las 
tecnicas de control quimico. 

En este trabajo se analiza por medio de un modelo matematico de sistemas de 
ecuaciones diferenciales Lotka-Volterra el efecto de la caza indiscriminada en el mar Adriatico 
en la poblacion de Selacios 1 (tiburones, mantarrayas, tintorera, etc.) durante la Primera Guerra 
Mundial 


1 “Se dice de los peces marinos cartilagineos que tienen cuerpo fusiforme o deprimido, cola heterocerca, pieI muy 
aspera, boca casi semicircular con numerosos dientes triangulares y de hordes cortantes o aserrados y mandibula 
inferior movil y varias hendiduras branquiales” (Real Academia Espanola). 
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ANTECEDENTES 


2.1. Ley Malthusiana 

Para conocer como evoluciona la poblacion de una determinada especie, incluyendo la 
humana, a lo largo de la historia se han considerado varios modelos. Uno de ellos, conocido 
como la Ley malthusiana 11 , que plantea el problema que representaba el creciente desequilibrio 
entre el incremental de la poblacion y el crecimiento de los satisfactores de esta. Malthus 
proponia que la razon de cambio de una poblacion era proporcional a la poblacion existente, es 
decir, si la variable x representaba la poblacion en cualquier instante t, donde t representa el 
tiempo, entonces se tiene 
dx 

— = kx donde k es la constante de proporcionalidad.(1) 

dt 

La ecuacion (1) es el primer modelo para el crecimiento de poblaciones; que 
corresponde a una ecuacion diferencial de variables separables y cuyo proceso de solution es: 

— = kdt — > f— = f kdt —» In x = kt + c ; 

x J x J 

que puede escribir como x(t) = ce kt .(2) 

Conviene senator que to ecuacion (2) solo se cumple para ciertas especies y por 
periodos determinados de tiempo. Asl, to Ley malthusiana es erronea para calcular el 
crecimiento de to poblacion de especies cuando el periodo de tiempo es grande, pues no 
considera factores tales como to extension del habitat, las epidemias, disponibilidad de 
alimentos, etc. Cuando se incorporan al modelo malthusiano limitaciones tales como el espacio 
para convivencia o to production de alimentos, solo por senator algunas variables, 
evidentemente to ecuacion diferencial cambia. 

Consideremos una vez mas to ecuacion propuesta por Malthus y consideremos que to 
constante k puede interpretarse como to diferencia entre los indices de natalidad y mortalidad 
a y b , respectivamente y que son funciones que dependen de to poblacion. Por lo que to 
ecuacion (1) puede reescribirse como 


n Propuesta por el economista ingles Thomas Robert Malthus en 1978 en su obra "An essay on the Principle of 
Population”. 
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(3) 


dx 

dt 


= (a-b)x 


2. 2. Ecuacion Logistica y funcion de Gompertz 

Si por razones o motivos tales como la extincion de especies, se considera que el Indice 
de natalidad se puede modelar con una funcion lineal decreciente, entonces se puede proponer: 
a = a 0 - a x x y sustituyendo en (3) tenemos: 

dx 

— = ( a 0 — a x x — b)x = (a 0 - b)x - a { x' 
dt 

Renombrando la diferencia de a 0 -b con una nueva variable tal que a 0 -b = K y haciendo 
a { =0,se obtiene otro modelo, llamado logistico 111 y que es muy usual para describir la razon 
de cambio, siendo esta 

— = Kx-^r 2 ; 
dt 

que puede tratarse como una ecuacion diferencial de variables separadas o como Bernoulli. 

En la decada de 1930, el biologo G. F. Gauze realizo un experimento con el protozoario 
Paramecium y uso una ecuacion logistica para modelar sus datos con resultados sorprendentes 
ya que la poblacion de protozoarios se comporta de manera muy similar a las predicciones del 
modelo. 

La ecuacion de Malthus y el modelo logistico no son las unicas ecuaciones propuestas 
para modelar el crecimiento de poblaciones, tambien existen modelos como la funcion del 
crecimiento de Gompertz y modelos de crecimiento dependiente de las estaciones como la 
ecuacion (4) para modelar cambios provocados por las estaciones en la disponibilidad del 
alimento, notese que las funciones periodicas son los elementos necesarios para modelar este 
tipo de comportamientos. 

— = kPcos(at-</>) .^ 

dt 


m La ecuacion diferencial logistica fue propuesta por el biologo matematico holandes Verhulst, como un modelo 
para el crecimiento de la poblacion mundial en la decada de 1840. 
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DESARROLLO 


3.1. La investigacion de D'Ancona 

En 1926 Humberto D’Ancona, un biologo italiano, completo un estudio estadistico de las 
poblaciones cambiantes de varias especies de peces de las orillas del norte del mar Adriatico. Su 
estimation de las poblaciones durante el periodo que va de 1910 a 1923 se baso en el numero 
de cada especie vendido en el puerto de Fiume, Italia. D’Ancona supuso, al igual que nosotros, 
que los numeros de varias especies en los mercados reflejaban la abundancia relativa de las 
especies en el Adriatico. En la tabla 1 se muestran los datos. 


Puerto/Ano 

1914 

1915 

1916 

1917 

1918 

1919 

1920 

1921 

1922 

1923 

Fiume 

12% 

21% 

22% 

21% 

36% 

27% 

16% 

16% 

15% 

11% 


Tabla 1 Porcentaje de selacios respecto al total de la pesca 


Como suele suceder, los datos no brindaron un apoyo abrumador para alguna teoria 
especifica de poblaciones cambiante de peces. Sin embargo, D’Ancona observo que los 
porcentajes de especies depredadoras fueron por lo general mayores durante e inmediatamente 
despues de la Primera Guerra Mundial. La pesca se redujo de forma considerable durante la 
guerra debido a que los Pescadores dejaron sus redes para pelear; D'Ancona concluyo que la 
reduction de la pesca dio lugar al cambio en las proporciones depredador-presa. Formulo la 
hipotesis de que durante la guerra la comunidad depredador-presa se aproximo a su estado 
natural de una proportion relativamente alta de peces depredadores, en tanto que la pesca mas 
intensa de los anos de la preguerra y la posguerra perturbaron el equilibrio a favor de las presas. 
Incapaz de explicar el fenomeno, D'Ancona pregunto a su suegro, el matematico Vito Volterra 
(1860-1940), si habia un modelo matematico que pudiera aclarar la cuestion. En unos meses 
Volterra describio una serie de modelos para las interacciones de dos o mas especies. 

3.2. Modelo depredador-presa (Lotka-Volterra) 

El modelo mas simple de la asociacion depredador-presa solo incluye el desarrollo o el 
deterioro natural y la interaction entre depredador y presa. Se supone que las demas relaciones 
son insignificantes. Volterra initio su analisis separando a la poblacion de selacios 
(depredadores), de los peces comestibles (presas): el poblacion de presas x(i) y la poblacion 
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de depredadores y(t) . Su razonamiento fue entonces que los peces comestibles no compiten 
muy intensamente entre si por su alimento, ya que este es muy abundante y la poblacion de 
peces no es muy densa. Por ello en ausencia de selacios, los peces comestibles crecerlan de 

dx 

acuerdo a la Ley Malthusiana —- = kx. Ademas la interaccion entre depredadores y presas se 

dt 

modela por los terminos de accion de masas proporcional al producto de las dos poblaciones 
axy slendo perjudicial para las presas; as! que la poblacion de peces comestibles en el tiempo 
t esta dada por: 
dx , 

— = kx- axy 
dt 

tambien concluyo que los depredadores tenian una tasa natural de decrecimiento - by, 
proporcional a su numero, y que su incremento tenia una tasa cxy .proporcional tambien a su 
numero, y al suministro de alimento x ; de esta forma la poblacion de selacios en el tiempo t 
esta dada por: 

dy , 

— = -by + cxy. 
dt 

De modo que el sistema de ecuaciones diferenciales El modelo para las poblaciones de 
depredadores y presas es el sistema depredador-presa (o tambien sistema llamado Lotka- 
Volterra lv ) 

dx dy 

— = kx- axy a — = -by + cxy .(5) 

dt dt 

describe la interaccion entre los selacios y los peces comestibles en el caso de no haber pesca 
alguna. A continuacion se estudiara este sistema y se obtendran algunas propiedades 
interesantes de sus soluciones. Despues, se incluira en el modelo el efecto de la pesca y se 
mostrara que un bajo nivel de la captura es mas benefico para los selacios que para las especies 
comestibles. De hecho, se Negara al sorprendente resultado de que un bajo nivel de pesca, en 
realidad, es danino para los peces comestibles. 


IV Volterra escribio “LesQons surla theorie mathematique de ta lute pour la vie” relacionado con sus teorlas en 1931. 
A. J. Lotka , un biologo estadounidense que despues fue actuario, llego de forma independiente a muchas de las 
conclusiones de Volterra; i /ease su libro, "Elements of Physical Biology”(1952). 


7 




3.3. Leyes de Volterra 


Volterra resumio sus conclusiones acerca de las soluciones y orbitas del sistema (5) en 
forma de tres leyes. Supongase que k,a,b,c son constantes positivas. 

3.3.1. Ley del Ciclo Periodico 

Ley del ciclo periodico. Las fluctuaciones de las poblaciones del depredador y su 
presa son periodicas. El periodo depende de los valores de los coeficientes de la 
tasa de cambio del sistema (5) y de los datos iniciales. Ademas, aumenta con la 
amplitud del ciclo correspondiente. 


Por lo tanto lo que tenemos que hacer es mostrar que el sistema (5) tiene una solucion 
periodica. Es decir, existe una constante T tal que x( t )= x( t + T ) y y( t )= (t + T) para toda t. 

A continuation se explica el proceso llevado a cabo para encontrar una ecuacion para las 
orbitas del sistema: 

Se divide la primera ecuacion de tasa de cambio del sistema (5) entre la segunda y se obtiene 

dx 

dt _ kx-axy _ kx axy _ kx axy 

dy -by + cxy -by + cxy -by + cxy y(—b + cx) y(—b + ex') 

dt 

De donde despues de simplificar el lado izquierdo y factorizar el termino derecho nos queda que: 

dx , x ..k 

— =(- - X—a) 

dy -by cx y 

Al separar las variables tenemos: 

(— — + CX )dx = (— — a)dy —> —-dx + cdx = —dy - ady 
x y x y 

Integrando ambos miembros de la ecuacion y pasando todos los terminos del lado derecho y 
multiplicando ambos miembros por (-1) tenemos: 

U I _ 7, f dy f J,, v (U 1„ „ , /7,1„ .. 


- b J— ycjdx = k J— -a^dy -> 
Donde C es una constante de integracion. 


{b In x - cx) y (k In y - ay) = C 


Elevando e a cada uno de los dos miembros 


b\n.x—cxyk\a.y—ay _ c 

b In x— cx k In y—ay s ^ 

- vy 


por la propiedad e = e e y proponiendo C= e c 


por la propiedad c In x = (In x) c obtenemos 
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(]nx) c -cx (]ny) k -ay _ s - In x _ 

v & — i-' por propiedad & — X 

(xV“)(/^) = C .(6) 

La formula (6) define una curva cerrada simple(es decir, un ciclo) para cada (X0,Y0) 
dentro del cuadrante de poblacion lo que significa que la solucion correspondiente 
x=x( t) & y= (y) del sistema es periodica.(vease llustracion 1). 



Ilustracion 1 Plano Fase 

De manera notable, la poblacion promedio de cada especie en cualquiera de los ciclos es 
una constante fija. Esta es la segunda ley de Volterra. 

3.3.2. Ley de los Promedios 

Ley de los promedios. En el sistema (5), las poblaciones promedio de depredador 
y presa durante un ciclo son, respectivamente, b/c y k/a 

Veamos por que se cumple la ley de los promedios. Supongase que x=x( t), y=y( t) es una 

solucion no constante que define un ciclo de periodo T. Las poblaciones promedio x y y 
durante un periodo se definen como 

1 T \ T 

x = — ]x(t)dt A y = — j \y(t)dt, 

^0 ^0 
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Demostraremos que x Primero se reordenan los terminos de la segunda ecuacion 
de tasa de cambio de (5) para obtener: 

y' / x _ b 1 y 

— = b + cx —> despejando X tenemos que: x v) ~ ~.(7) 

y c c y 

Luego integramos cada miembro de (7) de 0 a T 

X = -\x(t)dt = Ifeft- I r fI4ft = fc_l hy(D-l°yCO) = * 

Ti Tic T\c y - T 


Tic 


c T 


c 


0 0 

Lo anterior es posible debido a que y(T)=y(0). Con un argumento similar se demuestra que 


y = k I a y de esta manera queda demostrada la validez de la ley de los promedios. 


La ley de los promedios nos acerca mas a los datos reales antes mencionados para la 
captura de peces, ya que tales datos consisten en promedios. Sin embargo, aun no se toma en 
cuenta a los Pescadores. 

Por que la captura dana al depredador y ayuda a la presa 

La tercera ley de Volterra explica que sucede cuando son capturadas dos especies. El 
modelo mas simple es el de captura de esfuerzo constante, en el que la cantidad capturada por 
unidad de tiempo es proporcional a la poblacion. 

dx 

— = kx-axy - H^x-{k- H x - ay)x .(8) 

dt 

— = -by + cxy -H 2 y = (-b - H 2 + cx)y 
dt 

Los numeros negativos HI y H2 son los coeficientes de captura. Cuando hay captura, el 
punto de equilibrio que esta dentro del cuadrante de poblacion se desplaza a la izquierda y hacia 
arriba a partir de x=b/c, y=k/a al punto 

x = (b + H 2 )/c y = (k — H x )l a 

Se supone que H2< c. De otro modo, la captura intensa de la presa x no deja suficiente 
alimento para el depredador y, por consiguiente, esta especie tiende a la extincion. Por la ley de 
los promedios, los que corresponden a la poblacion alrededor de cualquier ciclo estan dados por 


10 





las coordenadas del punto de equilibrio. Puesto que la captura provoca que el punto de equilibrio 
se mueva hacia arriba y a la izquierda de la posicion original en el cuadrante de poblacion, la 
captura aumenta el promedio de la poblacion de presas pero disminuye el de la poblacion de 
depredadores. 

3.3.3. Ley de Captura 


Ley de captura. La captura de esfiierzo constante aumenta el numero promedio 
de presas por ciclo y disminuye el numero promedio de depredadores 


Antes de hacer comparaciones de los datos reales con las predicciones de la ley de 
Volterra es elemental reformular la ley en terminos de porcentajes( como los de la tablal) en vez 
de promedios. 


3.3.4. Ley de Porcentaje de Captura 


Ley de porcentaje de captura. La captura de esfiierzo constante aumenta el 
porcentaje promedio de las presas por ciclo en la poblacion total de peces y 
disminuye el porcentaje promedio de depredadores por ciclo. 


Si los coeficientes de captura en el sistema (8) son demasiado grandes, el punto de 
equilibrio interno (b + H 2 )/c,{k - A/j)/ acorta el eje x positivo y se extingue una o ambas 
especies. 

Con la ley de porcentaje de captura de Volterra se responde la pregunta de D'Ancona. Una 
disminucion en la tasa de captura da lugar a un incremento en el porcentaje de depredadores. 
(Vease la tabla 1). Cuando la tasa de captura disminuyo durante los anos de la guerra y luego 
aumento durante la posguerra, el porcentaje de depredadores aumento y luego disminuyo , justo 
como predice el modelo de Volterra. 
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Validez de la ley de captura de Volterra 

Desde su formulation, en muchas ocasiones el modelo de Volterra se ha puesto a 
prueba y recibido apoyo, pero continua siendo el punto de partida para muchos intentos serios 
de entender como evolucionan las comunidades de depredadores y presas con o sin captura. 

Una confirmation impresionante de la validez general de la ley de captura de Volterra 
ocurrio cuando se aplico el insecticida DDT V para controlar un hemiptero que infesto los huertos 
de cltricos en Estados Unidos de America. El insecto fue traIdo accidentalmente de Australia en 
1868. La cantidad de estos insectos fue controlada (pero no eliminada) mediante la importation 
de su depredador natural, un tipo particular de mariquita. Cuando se introdujo el DDT como 
agente de “captura”, se esperaba que el hemiptero pudiera ser eliminado por completo. Pero el 
DDT actuo de manera indiscriminada, matando el nurnero de mariquitas, en tanto que aumento 
la poblacion de hemipteros, liberados ya de la depredation por parte de las mariquitas. 


CONCLUSION 

La razon por la que durante la Guerra la tasa de poblacion de las presas disminuyo de 
manera radical se debe a que no sufrio de caza, generalmente en un entorno en el que el que 
se cazan presas, se pensarla que disminuye la poblacion de estas, pero dado que los 
depredadores dependen de las presas para aumentar su tasa de poblacion y que lo contrario, 
no sucede ,al haber pocos depredadores por ausencia de alimento las presas aumentan su tasa 
de poblacion exponencialmente, esto porque la tasa de natalidad de las presas es muy grande 
por ello la caza de presas dana al depredador y ayuda a la presa. 

Como lo expuso Volterra, podemos conjeturar que para un entorno natural (en el que 
la caza no es un factor a considerar) en ausencia de depredadores, la poblacion de las presas 
creceria exponencialmente tambien pero... ^Pueden subsistir en equilibrioambas especies? 

Obtuvimos que los puntos criticos del sistema son (0,0) y (b/c, k/a), los cuales fueron 
calculados a partir de la segunda ley de Volterra, y que ademas hacen a x’ =0 y y-0. El 
segundo conjunto de puntos criticos fue nuestro mayor interes. 


v El DDT(Dicloro Difenil Tricloroetano) fue utilizado con intensidad como insecticida, en los ultimos ahos se prohibio 
su danos a la salud humana 
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Si el punto de equilibrio del sistema tambien llamado punto(s) critico(s) queda dentro 
del cuadrante de poblacion, ademas, la curva que toca el punto (xo, yo) donde xo y yo son las 
poblaciones iniciales, este proxima al punto critico para un tiempo mayor a cero, importando que 
en ningun momenta las variables xo =0 o yo=0. Mientras mas alejado se encuentre un punto 
cualquiera (xo, yo), del punto de equilibrio no se puede asegurar con facilidad que ambas 
especies subsistan. 

Respecto al tiempo no podemos deducir las ecuaciones para x & y con respecto de t, pero 
podemos tener una relacion de la poblacion de y para cierta poblacion de x ,sabiendo que si en 
dado momenta para cierta poblacion de x se tienen 2 puntos de y es porque provienen de 2 
diferentes tiempos, en la realidad podriamos darnos cuenta de alguna manera cada cuanto se 
repite un periodo para poder establecer a x & y en funcion del tiempo. 
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